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Abstract 

The dynamics of the second order rational difference equation in the title with 
complex parameters and arbitrary complex initial conditions is investigated. Two 
associated difference equations are also studied. The solutions in the complex 
plane of such equations exhibit many rich and complicated asymptotic belabour. 
The analysis of the local stability of these three difference equations and period¬ 
icity have been carried out. We further exhibit several interesting characteristics 
of the solutions of this equation, using computations, which does not arise when 
we consider the same equation with positive real parameters and initial condi¬ 
tions. Many interesting observations led us to pose several open problems and 
conjectures of paramount importance regarding chaotic and higher order periodic 
solutions and global asymptotic convergence of such difference equations. It is our 
hope that these observations of these complex difference equations would certainly 
be new add-ons to the present art of research in rational difference equations in 
understanding the behaviour in the complex domain. 
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1 Introduction and Preliminaries 

Consider the difference equation 


a /9 

^n+l = I 

^n—1 


n = 0,1,... 


( 1 ) 


where the parameters a, (3 are complex numbers, and the initial conditions z_i and zq 
are arbitrary complex numbers. 


The same difference equation is studied when the parameters a and /3 and initial 
conditions are non-negative real numbers, Eq. 0 was investigated in [2] and [3j where the 
global asymptotic stability of the positive equilibrium was proved for all the parameters. 
In this present article it is an attempt to understand the same in the complex plane. 

The set of initial conditions z_i, Zq E C for which the solution of Eq. Q is well defined 
for all n > 0 is called the good set of initial conditions or the domain of definition. It is 
the compliment of the forbidden set of Eq. Q for which the solution is not well defined 
for some n > 0. See na for the definition and further work on obtaining the forbidden 
set for the first order Riccati Difference equation. However, for the second and higher 
order rational difference equations, the lack of an explicit form for the solutions, makes it 
very challenging to obtain the good set. For Eq. 0 we pose the following open problem 
of paramount importance and difficulty. 

Open Problem 1.1. Determine the good set of initial eonditions for Eq. 0 

Our goal is to investigate the character of the solutions of Eq. 0 when the parameters 
are real or complex and the initial conditions are arbitrary complex numbers in the 
domain of definition. For the rest of the sequel we assume that the initial conditions are 
from the good set na. 

We now present some preliminary material which will be useful in our investigation 
ofEq.0. 

Let / : —>■ C be a continuous function, where /C is a non-negative integer and 

C is an interval of complex numbers. Consider the difference equation 

^n+l f ^n—11 • • • t ^u—k) ) ^ 0 , 1 ,... (2) 

with initial conditions Z-ic, Z-ic+i, ..., zq E C. 

We say that y is an equilibrium point of Eq. Q if 

/(z,z, ...,z) = z. 

We now impose the further restriction that the function f{uo,ui,...,uic) be 
continuously differentiable. 
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The linearized equation of Eq.([^ about the equilibrium ^ is the linear difference 
equation 


^n+l T ^l^n—1 T ' ' ' T ^K^n—K i ^ 0, 1, . . . (3) 

where for each i = 0,1,..., /C 


df 

ai = 

The characteristic equation of Eq. Q is the eqnation 

— aoA^ — aiA^“^ _ ... _ ax;-iA — a^c = 0. (4) 

The following result, called the Linearized Stability Theorem, is usefnl in determining 
the local stability character of the equilibrium A of Eq. (§. n. 

Theorem A (The Linearized Stability Theorem) 

The following statements are true: 

1. If every root of Eq.^ has modulus less than one, then the equilibrium z of Eq.^ 
is locally asymptotically stable. 

2. If at least one of the roots of Eg. Q) has modulus greater than one, then the 
equilibrium z of Eq. ^ is unstable. 

The equilibrium solution A of Eq.Q is called hyperbolic if no root of Eq.Q has 
modulus equal to one. If there exists a root of Eq. Q with modulus equal to one, then 
A is called non-hyperbolic. 

The eqnilibrium point A of Eq.([^ is called a sink if every root of Eq.Q has modnlns 
greater than one. 

The equilibrium point f of Eq.Q is called a saddle point equilibrium if it is hyper¬ 
bolic, and if in addition, there exists a root of Eq.Q with modulns less than one and 
another root of Eq.Q with modulus greater than one. In particular, if A is a saddle 
point equilibrium of Eq.([^, then A is unstable. 

The equilibrium point A of Eq.([^ is called a repeller if every root of Eq.Q has 
absolute value less than one. 

The following theorems would be useful in determining the characteristics of zeros 
of the characteristic polynomial. 
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Theorem B (Enestrom-Kakeya Theorem) 

Consider a polynomial P{z) = Yl]=o degree n with real coefficients ao, ai,... 

such that ttn > ttn-i > a„_i • • • > Uq > 0, then P{z) has all its zeros in for | 2 |< 1 
[8]. In literature [9], there are some extensions and generalizations of this theorem. In¬ 
terestingly, in 1967, Govil and Rahaman [5] extended it to a polynomial with complex 
coefficients by proving the foflowing theorem. 

Theorem C (Govil-Rahaman Theorem) 

Let P{z) = ^ polynomial of degree n with complex coefficients oq, oi,... a„, 

there exists an a > 0 such that 

\ ttn \> a \ 0„_1 |> I ttn-l \> ■ ■ ■ > 0 ”“^ I Q-1 |> a"" I Go I 

then the polynomial P{z) has all its zeros in | z |< ^Ri where Ri is the greatest positive 
root of the cubic trinomial 7?’*+^ — 27?” -|- 1 = 0. 


2 Local Stability of the Equilibriums 


In this section we establish the local stability character of the equilibria of Eq. Q when 
the parameters a and (5 are considered to be complex numbers with the initial conditions 
are arbitrary complex numbers. 

The equilibrium points of Eq. Q are the solutions of the equation 

_ a B 

z = — — 

z z 

Eq.Q has two equilibria points 

Zi,2 = 

The linearized equation of Eq.Q with respect to the equilibrium zi 2 is 

a „ jB 


Zn+l + 


-,Zn + 


01 (B oi (B 

with associated characteristic equation 


Zn—l 0 , fl 0 , 1 ,..., 


A2 + 


a 


;A + 




= 0 . 


jB Oi jB 

The following result gives the local asymptotic stability of the equilibrium Zi^ 2 - 


( 5 ) 
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Lemma 2.1. The equilibriums zi ^2 = ±-\/a; + /5 of -Eg. Q) are locally asymptotically 

i ^ for some a > 0, f denotes golden ratio. 


y 

< 

a 

ol^(3 


a+/3 


stable if and only if a 


Proof. From Theorem C, the all roots of the characteristic equation (5) he in the close 
ball \z\ < \Ki where Ki is the largest zeros of the cubic polynomial — 2K^ + 1 such 
that for some a > 0 


1 > o 


a 

> 

3 

OL ^ 

OL P 



( 6 ) 


is satisfied. 

In this case 


Ki = 


\/5 + l 


which is known also as golden ratio, f. Therefore the closed ball becomes 

z\ < — 


(7) 


From Theorem (A) we have if every roots of the Eq. (5) has modulus less than one, 
then the equilibriums Zi ^2 are locally asymptotically stable. Then from Eq. (7), the radius 
of the ball \z\ < - should be - < 1. In turns, it is o > 6. 

II — n n ' I 


Therefore, The equilibriums of the Ec 


asymptotically stable if and only if a 




a+y 


■ (6) Zi ^2 = ±-\/« + Id of Eq. § are locally 


< 


ct+/3 




□ 


Lemma 2.2. The equilibriums zi ^2 — ±-\/a; + (d of Eg. Q) are sink if and only if 
< - for some 0 < a < f denotes golden ratio. 


y 

< 

a 

a+/3 


a+(3 


a 

Proof. From the Theorem (A), we have if every roots of the Eq. (5) has modulus 
greater than one, then the equilibriums Zi ^2 are sink. Then from Eq. (7), the radius of 
the ball \z\ < - should be - > 1. In turns, it is a < cb. 

Therefore, The equilibriums Zi ^2 = 2sz^a + j3 of Eq.(l) are sink if and only if a EE < 
< - for some 0 < o < </>. 


0.^(3 


a+/3 


□ 


Lemma 2.3. The equilibriums Zi ^2 = ±\/« + 3 of Eq.^ are hyperbolic if and only if 

— Q: j ( —OL \2 A /3 ~ 

2 


^1. 


Proof. The two roots X± (say) of the Eq. (6) are 

— OL 

ct+/3 


A±(q;,/5) = 


± - 4 


_±_ 

a+y 
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From the Theorem (A), the equilibriums zi ^2 are hyperbolic if the modulus of the roots 
of the Eq. (6) are not equal to 1. If there exists a root of Eq. (6) with modulus equal 
to IS, then the equilibrium is called non-hyperbolic. So, if the modulus of X± should be 
non equal to 1 then the equilibriums Zi ^2 are hyperbolic. That is, |A±| 1. 

In addition to the hyperbolicity of the equilibriums, if |A+| > 1 and |A_| < 1 then the 
equilibriums Zi ^2 are saddle. 

If any one of the parameters a or /3 = 0, then |A+| = 1 then the equilibriums Zi ^2 are 
non-hyperbolic. 

□ 

In different specific cases of the parameters a and /5 of the Eq. (1), the stability of 
the equilibriums \/a-\- jd are characterized through the characteristic equation Eq. (5) 
about the equilibriums of the Eq. (1) are described in the Table [1], 


Parameters 

Characteristic Equation 

Modulus of the Zeros 

Inference 

a = (3 

A^ + ^A+^=0 

1 A± |> 1 

±\/^ are sinks 

a = il3 

A^ + T+? = 0 

1 A± |< 1 

±\//3 + are La.s 

a = —if3 

A = o 

1 A± |> 1 

±\/^ are sinks 

a = 0 & d ^ 0 

A2 + 1 = 0 

1 A± h 1 

±v^ are non-hyp 

Oi 0 & /3 = 0 

A2 + A = 0 

1 A± hO,l 

3z^/a are non-hyp 


Table 1: Local stability of the equilibriums of Eq. (1) 


3 Dynamics of Associated Difference Equations 

Here we consider another two associated difference equations which are defined as follows 


1 R 

— Cf + P 

^n—l 

n = 0,1,... 

(8) 

^n—1 . n 

Zn+i = a + /5 

n = 0,1,... 

(9) 
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where the parameters a, (3 are complex numbers, and the initial conditions Z-i and zq 
are arbitrary complex numbers. 

These two difference equations Eq. (8) and Eq. (9) are derived from Eq. (1) just by 
multiplication of Zn and Zn^i with the function argument respectively . 

It is natural to pose a similar open problem as stated before for the Eq. (1). 

Open Problem 3.1. Determine the good set of initial eonditions for Eq.(8) and Eq. 

(9). 

Let us now investigate the local stability of the equilibriums of these two difference 
equation in the following section. 


3.1 Local Stability of the Equilibriums 

In this section we describe the local stability character of the equilibria of Eq.(8) and 
Eq.(9) when the parameters a and (3 are considered to be complex numbers with the 
initial conditions are arbitrary complex numbers. 

The equilibrium are for the difference equations Eq.(8) and Eq.(9) is a -\- (3. Let us 
linearize the Eq.(8) and Eq.(9) about the equilibrium « + /5 as follows. The characteristic 
equations associated to the linear equations corresponding to Eq.(8) and Eq.(9) are 


A' 


tt + p tt + p 


0 . 


( 10 ) 




a 


:A- 


a 


(X jS (X jS 


= 0 . 


( 11 ) 


In different specific cases of the parameters a and (3 of the Eq. (8), the stability of the 
equilibriums a-\- (3 are characterized through the characteristic equation Eq. (10) about 
the equilibriums of the Eq. (8) are described in the Table [2]. 
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Parameters 

Characteristic Equation 

Modulus of the Zeros 

Inference 

a = (3 

f A + 1 = 0 

1 A± |< 1 

2a is l.a.s 

a = ij3 


1 A± |> 1 

P + ip is sinks 

a = —ip 

a2-t^a + ^ = o 

A± > 1 

P — ip is sinks 

oi = 0 & /3 ^ 0 

A2 - A + 1 = 0 

1 A± h 1 

P are non-hyp 

Qi ^ 0 & /3 = 0 

A2 = 0 

1 A± 1= 0 

a are repeller 


Table 2: Local stability of the equilibriums of Eq. (8) 

In different specific cases of the parameters a and /5 of the Eq. (9), the stability of the 
equilibriums a + /5 are characterized through the characteristic equation Eq. (11) about 
the equilibriums of the Eq. (9) are described in the Table [3]. 


Parameters 

Characteristic Equation 

Modulus of the Zeros 

Inference 

a = P 

A2+iA-i=0 

\ ^ 

2a is non-hyp 

a — ip 

- lil = 0 

1 A± |< 1 

P -\- ip is l.a.s 

a = —ip 


1 A± h 0 

P — ip is repeller 

Qi = 0 & /3 ^ 0 

A2 = 0 

T 

0 

P are repeller 

(a ^ 0 & /3 = 0 

A2 + A - 1 = 0 

1 A± |> 1 

a are sinks 


Table 3: Local stability of the equilibriums of Eq. (9) 



4 Periodic of Solutions 


In this section we discuss the global periodicity and the existence of solutions that con¬ 
verge to periodic solutions of Eq.(l), Eq.(8) and Eq.(9). 

A difference equation is said to be globally periodic of period t if for any 

given initial condition. The characterization of periodic rational difference equations (ei¬ 
ther with real or complex coefficients) of order k is a challenging area of current research. 
See [1] and [H] for several periodicity results of second and third order rational difference 
equations, respectively, with nonnegative real coefficients and arbitrary nonnegative real 
initial conditions. 


4.1 Periodic Solutions of period 2 and 4 

Eor the difference equation Eq.(l), it is our firm conviction that for a = 0 or /3 = 0, the 
all solutions are periodic and of period 4 or 2 respectively. Also, it is conjectured that 
all solutions of the difference equation Eq.(l) with the parameter /5 = 0 or a = 0 are 
convergent and convergent to a periodic point or of period 2 or 4 respectively. 

Eurther we took an attempt to investigate the higher order periodic cycles. We found 
there there is no periodic solutions of period 3 which is formally proved in the next 
section. Interestingly there are periodic solutions of higher periods viz. 2,4,5,6,7,.... 


4.2 Higher Order Cycles in Solutions 

In investigating the solutions of higher order cycles of order d of the Eq. (1) the following 
system of nonlinear equations needs to be solved. The solutions are indeed the cycles. 

oi jd , , 

^mod{n,d) T 

^mod{n—l,d) ^mod{n—2,d) 

for n = 1,2,3,..., d and where mod(n, d) denotes the x such that n = x mod{d). 


Theorem D: There does not exist any periodic solution of period 3. 
Proof: The solutions of the equations 


f a P a P a 

SZq— - 1 - ,Zi = - 1 - ,^ 2 =- 1 ” 

( 2:2 Zi Zo Z 2 Zi 


I 

Zo 


are the periodic solutions of period 3. Interestingly there is no solutions except the 
equilibrium of the equation Eq. (1). Hence there is no solution of period 3. 


A list of difference periodic cycles of different length are adumbrated for the param¬ 
eters a = P = 1 m the Table 4. 
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Parameters 

CL 

One of the Solutions 

a = l;/3 = 1 

4 

0.765367h -1.84776h -0.765367^ 1.84776f 

Of = 1; /3 = 1 

5 

0.309721i, -1.83083i, -2.68251i, 0.918986i, -0.71537i 

a — 1] (3 = 1 

6 

0.53713i, -0.735107i, -0.501402i, 3.35475i, 1.69632i, -0.887595i 


7 

0.563218h -1.45984f, -1.0905H, 1.6020H, 0.29279R, -4.03963f, -3.16786f 


Table 4: Cycle solutions of different order for different choice of a and /5 of Eq. (1), 
Here CL denotes cycle length. 


There are also many higher order cycles too. For a = (15, —88) and P = (—53, —30) and 
the initial values Zq = (65,-97) and Zi = (-92,-67), we have encountered a solution 
of the difference equation Eq. (1) which converges to a periodic point of period 23. 
The periodic cycles is (18.574,-4.5796) ^ (3.295,-1.1914) ^ (-9.4511,-17.474) ^ 

(15.164,-23.792) ^ (3.424,0.61574) ^ (-13.604,-7.0189) ^ (3.742,-25.505) ^ 

(2.618,3.3548) ^ (-9.7582,5.5619) ^ (-11.358,-10.844) ^ (-1.2796,5.1965) ^ 
(0.10378,15.669) ^ (-18.574,4.5796) ^ (-3.295,1.1914) ^ (9.4511,17.474) ^ 

(-15.164,23.792) ^ (-3.424,-0.61574) ^ (13.604,7.0189) ^ (-3.742,25.505) ^ 

(-2.618,-3.3548) ^ (9.7582,-5.5619) ^ (11.358,10.844) ^ (1.2796,-5.1965) ^ 
(—0.10378, —15.669) —>■ (18.574, —4.5796). The plot of the real and imaginary sequence 
is given the Fig. (1). 


Plot of Real sequence 



100 200 300 400 500 600 700 800 

Iterations 

Plot of Imaginary sequence 



Figure 1: Plot of Real and Imaginary Sequence for a = (15, —88) and /3 = (—53, —30) 
and the initial values Zq = (65, —97) and Zi = (—92, —67). 


We also have listed different periodic cycles of the different length for the associated 
difference equations Eq. (8) and Eq. (9) in the Table 5 and Table 6. 
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Parameters 

CL 

One of the Solutions 

a = l;/3 = 1 

3 

1.24698, -1.80194, -0.445042 

Of = 1; /3 = 1 

5 

0.83083, -1.91899, -1.30972,1.68251, -0.28463 

Of = 1; /3 = 1 

6 

1.80194, -0.445042,1.24698, -1.80194, -0.445042,1.24698 

a = l;/3 = 1 

7 

1.87278, -0.556474,0.702862, -0.263063,0.625725, -1.37861, -1.20322 


Table 5: Cycle solutions of different order for different choice of a and /5 of Eq. (8). 


Parameters 

CL 

One of the Solutions 

Of = 1; /3 = 1 

3 

1.24698, -0.445042, -1.80194 

a = l;B = l 

4 

1.53339 - 0.608009i, 1.81536 + 0.929423b 1.53339 - 0.608009b 1.81536 + 
0.929423f 

Of = 1; /3 = 1 

5 

0.574313 -6 0.798528i,-0.273032 - 0.160806i,-1.84063 - 

1.25163b 1.14206 - 0.00923437i, -0.60271 - 1.1089i 

a = l;/3 = 1 

6 

2.61506,1.61917,2.61506,1.61917,2.61506,1.61917 

Of = 1; /3 = 1 

7 

0.962688 + 0.453798b-0.251383 - 0.177869i,-2.40312 + 

0.602711i, 1.08095 0.0943187i,-1.15807 0.745878i, 0.377343 - 

0.482479b -1.12397 - 0.739107i 

a = l;/3 = 1 

8 

27.0466,1.03839,27.0466,1.03839,27.0466,1.03839,27.0466,1.03839 


Table 6: Cycle solutions of different order for different choice of a and P of Eq. (8). 


Theorem E: There does not exist any periodic solution of period 4. 
Proof: The solutions of the equations 


^1 


, Pzo 

Ot H-, Z 2 

Z3 


, Pzi PZ2 , PZ3\ 

a. H- , Zs — O' -\ - ^ zq — O' -\ -> 

Zo Zi Z 2 J 


are the periodic solutions of period 4. There is no solutions except the equilibrium a + /f 
of the equation Eq. (8). Hence there is no solution of period 4. 


Theorem F: The 2-cycle of the difference equation Eq. (1) are attracting, repelling and 
non-hyperbolic if and only if the modulus of the eigenvalues of the matrix 

a /3 

(3—a 13—a 

— a[3 
(a-yp 


are lesser, greater than or equal to 1 respectively. 


Proof: The two cycles of the difference equation Eq. (1) are the solutions of the nonlinear 
equations 

a B a B 

Zi = - 1 - , Zq = - 1 - 

Zq Zi Zi Zq 
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Therefore the 2 cycles are —\/a — (3 and ^/a + /3 which is actually hxed point of the Eq. 
(1). So we are mainly interested in the other 2 cycle which is —y/a — jS. 


The characteristic equation of the linearized equation about the 2-cycle —\/a — {5 of 
Eq. (1) is 


ct (3 

j3—OL (3—a 

— q :^+/ 3 q ;+/ 3 ^ ocj3 

(q:-/9)2 {a-PY 


= 0 


matrix 


The zeros of this characteristic equation, which are essentially the eigenvalues of the 

[3—a 


(3—0. 

—q;^+/3q;+/3^ 

{(y--l3Y 


a(3 


are 


—4q;^/5 — 8q;^/5^ -|- 4/5^ -|- 2q;/5 —o? -|- \J— 8q;^/5^ -|- 4/5^ -|- 2q;/5 


2{a — py 


2{a-yy 


The modulus of these two eigenvalues need to be lesser, greater or equal to 1 in de¬ 
termining the 2-cycle as attracting, repelling or non-hyperbolic. Hence the theorem is 
proved. 


In the similar fashion, the local stability of the 2-cycle of the difference equations Eq. 
(8) and Eq. (9) can also be analyzed easily which we left to reader. This result would 
encourage us to draw the following conjecture. 

Conjecture 4.1. For every finite length (if exists) cycle of the difference equations Eq. 
(1), Eq. (8) and Eq. (9), there exists at least one solution which converges to the cycle. 


5 Chaotic Solutions 

For the difference equations Eq. (1), Eq. (8) and Eq. (9) we have found chaotic solu¬ 
tions. The chaotic property of the solutions can be ensured through the largest positive 
Lyapunav exponent of the solutions. Few examples are illustrated in the following table 
(7) with corresponding plots of the solutions as shown in Fig. (2). The green and blue 
plot are denoting real and imaginary sequences respectively. 


Zn+1 — ^ ^ 

Zn+i - a + 13 

Zn+l - 

a = (30,47), 13 = (30,-10), 
zq = (9, -41), zi = (49. - 63) 

a = (56, -22), /? = (-52, -19), 
zo = (-81,-74), zi = (89,92) 

a = (4,-81), /? = (64,64), 
zo = (45,-70), = (32,4) 


Table 7: Cycle solutions of different order for different choice of a and fi of Eq. (8) and 
initial values. 
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Figure 2: Chaotic Solutions for the difference equations Eq. (1) (Top), Eq. (8) (Middle) 
and Eq. (9) (Bottom). 


The Lyapunav exponent of the solutions of the three cases as stated in the table are 
1.6015, 1.2414 and 0.6885 respectively which declare that the solutions are indeed 
chaotic. 

Open Problem 5.1. Find out the set of all possible initial values Zq and Zi for which 
the solution are chaotic for given parameters a and jd. 

6 A Comparative View of Dynamics 

In this section, an attempt has been made to compare the dynamics of the three asso¬ 
ciated differences equations Eq. (1), Eq (8) and Eq (9). The comparison of the limiting 
behaviour of the dynamics (sequence of iterates) and the state spaces over 50000 iter¬ 
ations of the three difference equations are presented in the following two subsections 
through examples. 

6.1 Comparison of Characteristics of Dynamics 

Here we present a table of comparative dynamics among the three difference equations 
Eq. (1), Eq (8) and Eq (9). It turns out that all these three difference equations are 
mostly having different characteristics for same parameters a and (d and with the initial 
values zq and zi. The fixed point of the difference equation Eq. (8) and Eq. (9) is 
a -\- jd. We observed that given the fixed parameters a and Id and the intimal values, 
any solutions of Eq. (8) which is convergent and converges to the the fixed point a + Id 
fetch another solution which is also convergent and converges to the same fixed point of 
the other difference equation Eq. (9). 
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Initial Values 

Zn^l = — H- — 

Zn Zn-1 

^n+1 —0^+ B 

Zn+l — 

a = (30,47), 

/I = (30,-10), 
zo = (9,-41), 
zi = (49, -63) 

Chaotic 

Convergent and con¬ 
verges to (60,37) 

Convergent and con¬ 
verges to (60,37) 

a = (56,-22), 

/3 = (-52,-19), 

ZQ = 

(-81,-74), 
zi = (89,92) 

Converges to a 

periodic point 

(10.393,-0.14432) 
of period 2 

Chaotic 

Converges to a periodic 
point ((33.8,60.46) of 
period 2 

a = (4,-81), 

/3 = (64,64), 
zo = (45,-70), 

= (32,4) 

Converges to a periodic 
point (6.9614,-10.414) 
of period 2 

Divergent 

Chaotic 

a = (15, —88), 

P = (-53, -30), 
zo = (65, -97), 
zi = (-92,-67) 

Periodic and of Period 

23 

Convergent and con¬ 
verges to (—38, —118) 

Convergent and con¬ 
verges to (—38, —118) 


Table 8: Local stability of the equilibriums of Eq. (9) 


The largest Lyapunav exponent of the solutions of the difference equations Eq. (1), Eq. 
(8) and Eq. (9) of the first three cases have already calculated in the previous section 
which ensure that the solutions are chaotic. 

6.2 Comparison of State Spaces of Dynamics 

In this section we have adumbrated a list of cases to compare state spaces over 50000 
iterations of the difference equations Eq. (1), Eq. (8) and Eq. (9) in the Table 9. It 
turns out that no state space is similar to others. 
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Initial Values 

Zn^l = — H- — 

Zn Zn-1 

—0^+ B 

Zn — l 

Zn+l — 

a — (9, —73), 

/3 = (-70, -49), 
zo = (52,110), 

zi = (68,88) 

Convergent and 

converges to 

(10.591,-5.759) 

Unbounded 

Convergent and con¬ 
verges to (79, —122) 

a = (100,-55), 

/3 = (31,21), 

zo = (-82,160), 

= (-11,-94) 

Unbounded 

Convergent and con¬ 
verges to (131, —34) 

Convergent and con¬ 
verges to (131, —34) 

a = (-29,33), 

/3 = (-44, -54), 

ZQ = (152,122), 

= (87,-191) 

Fractal 

Unbounded 

Convergent and con¬ 
verges to (—73, —33) 

a = (58,56), 

/3 = (34,-74), 
zo = (-8, -59), 
zi = (-57,-91) 

Unbounded 

Unbounded 

Convergent and con¬ 
verges to a periodic 
point (56.376,-118.56) 

a — (98,1), 

/3 = (-46, -80), 
zo = (-99,130), 

= (55,75) 

Fractal 

Convergent and con¬ 
verges to a periodic 
point (117.82,14.575) of 
period 6 

Convergent and 

converges to a 

periodic point 

(-61.076,-143.49) 
of period 2 

(X = (—64,0), 

/I = (4,99), 
zo = (-89,184), 

= (29, -32) 

Unbounded 

Unbounded 

Convergent and con¬ 
verges to (—60, 99) 

a = (80, —87), 

(3 = 

(-33,-100), 
zo = (87,64), 

^1 = (42,49) 

Fractal 

Chaotic 

Convergent and con¬ 
verges to (47, —187) 

a = (4,55), 

/3 = (-76,25), 
zo = 

(-34,-32), 
zi = (64,6) 

Unbounded 

Unbounded 

Convergent and con¬ 
verges to a periodic 
point (19.794,105.86) of 
period 2 


Table 9: Characterizations of State space over 50000 iterations of Eq. (1), Eq. (8) and 
Eq. (9) 


The states spaces are hgured corresponding to each of the above six cases in the Fig. 
(3). The states space is the set of all Zi for all i ranges from 1 to 50000. 
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Figure 3: States space of 50000 iterations of the Difference equations Eq. (1), Eq. (8) 
and Eq. (9). Each row of the figures denotes each different case as states in the Table 
(9) respectively. 
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In the case of 3rd, 5th and 7th cases of the Table (9), the solntions of the difference 
equation Eq. (1) are turned out to be fractal. The complex sequence plot of the different 
number of iterations are given in the in Fig. (4) for the 3rd and 5th cases. 


Plo of Complex sequence 



Plot of Complex sequence 







Plot of Complex sequence 



Plot of complex sequnce 







1.6 1.6 2 


Plot of complex sequence 



Plot of Complex sequence 




Iterations 

Plot of Complex sequence 




Iterations 



Iterations 


Figure 4: Plot of Complex sequences which turns out to be a fractal over different 
number of iterations; First two rows and next consecutive two rows of figures for 3rd 
and 5th cases respectively 
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The fractal dimension of the fractal solution for the 3rd, 5th and 7th cases are 1.82779, 
1.89333 and 1.9127 respectively. This fractal dimension have been calculated through 
the software Benoit^^. 

7 Conclusion 

Nonlinear difference equations are an established dynamic area of research in both real 
and complex domain, but the study of rational difference equations is still in its infancy. 
The first order Riccati difference equation, for which the solution could be obtained 
explicitly, has been studied both in the real and complex domains, see ini and na. 
Second and third order rational difference equations with positive real parameters and 
initial conditions, have also been studied extensively, see [Ij. Also see [Hj and ng. 

When we consider the parameters and initial conditions to be complex, there is 
hardly anything known about rational difference equations. In the present manuscript, 
we took three closely related difference equations and did lot of computational study in 
understanding their dynamics characteristics. It is turned out as expected the richness 
and complexity of the dynamics in terms on chaotic, fractal, unbounded solutions where 
in the case of real parameters and real initial values nothing were found for the same 
difference equations. Our main purpose behind this manuscript is to probe and analyze 
how the dynamics changes for rational equations when we shift from positive real domain 
to the complex domain. 

We have carried out some local stability and global periodicity analysis. We have 
observed computationally the different characteristics of boundedness, unboundedness, 
convergence, convergence to periodic solutions, fractal and chaos. Each one of these 
characteristics of the dynamics is important in their own right and we have posed them 
as challenging open problems which require further theoretical investigations. Further, 
certain peculiar behavior of the dynamics in the complex domain has been brought 
forward which is in stark contrast to the positive real scenario. 

The work carried out is generic in nature and can be helpful towards the study of 
other rational difference equations of second and higher orders. It is our believe, that 
this work would initiate the study and understanding of rational difference equations in 
the complex plane. 
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